Frebenius method